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We investigate the effects of bosonic trilinear interactions in preheating after chaotic inflation. A 
trilinear interaction term allows for the complete decay of the massive inflaton particles, which is 
necessary for the transition to radiation domination. We found that typically the trilinear term is 
subdominant during early stages of preheating, but it actually amplifies parametric resonance driven 
by the four-legs interaction. In cases where the trilinear term does dominate during preheating, the 
process occurs through periodic tachyonic amplifications with resonance effects, which is so effective 
that preheating completes within a few inflaton oscillations. We develop an analytic theory of 
this process, which we call tachyonic resonance. We also study numerically the influence of trilinear 
interactions on the dynamics after preheating. The trilinear term eventually comes to dominate after 
preheating, leading to faster rescattering and thermalization than could occur without it. Finally, 
we investigate the role of non-renormalizable interaction terms during preheating. We find that if 
<**] ' they are present they generally dominate (while still in a controllable regime) in chaotic inflation 

models. Preheating due to these terms proceeds through a modified form of tachyonic resonance. 



I. INTRODUCTION 



According to the inflationary scenario, the universe at early times expands quasi-exponentially in a vacuum-like 
' state without entropy or particles. During this stage of inflation, all energy is contained in a classical slowly moving 
t^J- , inflaton field. Eventually the inflaton field decays and transfers all of its energy to relativistic particles, which starts the 
i-H ' thermal history of the hot Friedmann universe. Particle creation in preheating, described by quantum field theory, 
is a spectacular process where all the particles of the universe are created from the classical inflaton. In chaotic 
inflationary models, soon after the end of inflation the almost homogeneous inflaton field coherently oscillates with a 
very large amplitude of the order of the Planck mass around the minimum of its potential. Due to its interactions 
with other fields, the inflaton decays and transfers all of its energy to relativistic particles. If the creation of particles 
,-SZh ' is sufficiently slow (for instance, if the inflaton is coupled only gravitationally to the matter fields) the decay products 
simultaneously interact with each other and come to a state of thermal equilibrium at the reheating temperature 
T r . This gradual reheating can be treated with the perturbative theory of particle creation and thermalization 
However, for a wide range of couplings the particle production from a coherently oscillating inflaton occurs in the non- 
1 perturbative regime of parametric excitation Q,y, 4, 5j. This picture, with variation in its details, is extended to other 
inflationary models. For instance, in hybrid inflation (including D-term inflation) the inflaton decay proceeds via a 
tachyonic instability of the inhomogeneous modes which accompany the symmetry breaking Q . One consistent feature 
of preheating - non-perturbative copious particle production immediately after inflation - is that the process occurs 
far away from thermal equilibrium. The transition from this stage to thermal equilibrium occurs in a few distinct 
stages, each much longer than the previous one. First there is the rapid preheating phase, followed by the onset of 
turbulent interactions between the different modes. Our understanding of this stage comes from lattice numerical 
simulations 7] as well as from different theoretical techniques @- For a wide range of models, the dynamics of scalar 
field turbulence is largely independent of the details of inflation and preheating 9]. Finally, there is thermalization, 
ending with equilibrium. In general, the equation of state of the universe is that of matter when it is dominated by the 
coherent oscillations of the inflaton field, but changes when the inflaton decays into radiation-dominated plasma |10| . 

Most studies of preheating have focused on the models with <f> 2 x 2 four-legs interactions of the inflaton <f> with 
another scalar field \- A common feature of preheating is the production of a large number of inflaton quanta with 
non-zero momentum from rescattering, alongside with inflatons at rest. The momenta of these relic massive inflaton 
particles eventually would redshift out. However, the decay of inflaton particles through four-legs <jxj) — > XX processes 
in an expanding universe is never complete. Thus inflaton particles later on will have a matter equation of state and 
come to dominate the energy density, which is not an acceptable scenario. Therefore, to avoid this, we must include 
in the theory of reheating interactions of the type (j>x n > that allow the inflaton to decay completely, thus resulting in 
a radiation dominated stage. Trilinear interactions are the most immediate and natural interactions of this sort. 
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Trilinear interactions occur commonly in many theoretical models. Yukawa couplings, for example, lead to trilinear 
vertices with fermions. Three-legs decay via fermions was in fact the first channel of inflaton decay considered in 
early papers on inflation pj within perturbation theory. It was later realized that even in the case of interaction 
with fermions the inflaton decay typically occurs via non-perturbative parametric excitations fill ]. However, three- 
legs decay via infractions with bosons is expected to be a dominant channel (due to Pauli blocking for fermions). 
Gauge interactions lead to trilinear vertices with vector fields. Even if we restrict ourselves to scalar field interactions, 
trilinear interactions naturally appear in many contexts. 

Consider for instance a chaotic inflation model with the effective potential V — zt^m 2 <j) 2 + j\(f) A + \ (j> 2 x 2 ■ The 
negative sign corresponds to spontaneous symmetry breaking <j> — ► <j> + a , which results in a classical scalar field VEV 
a = -^7=. The interaction term g 2 (ft 2 x 2 then gives rise to a trilinear vertex ag 2 4>\ 2 along with the four- legs interaction. 
Spontaneous symmetry breaking naturally emerges in the new inflationary scenario. However, this results in a mass 
for the x particles which is comparable to the inflaton mass. This complicates the bosonic decay of the inflaton after 
new inflation, see also |12| . One encounters the same problem for spontaneous symmetry breaking in chaotic inflation 
models. However, a very different picture appears for chaotic inflation in supersymmetric theories. 

Trilinear interactions occur also naturally in supersymmetric theories. Typical superpotentials are constructed from 
trilinear combinations of superfields, W = J2ijk ^ijk&i&j&k- The cross terms in the potential V — J2i \Wi\ 2 then 
give rise to a three-legs interaction, in addition to a four-legs vertex. Besides many other advantages, supersymmetry 
is useful to protect the flatness of the inflaton potential from large radiative corrections. Therefore we will pay special 
attention to trilinear interactions in SUSY theories. In SUSY the strengths of the couplings of three-legs and four-legs 
interactions are tightly connected. 

The scalar field potentials that we consider in this paper, and their stability constraints, are discussed in Section UTI 

In general, we expect the inflaton to decay simultaneously via all potential interactions: three-legs, four-legs etc. In 
many cases, notably in the supersymmetric case, the leading channel of inflaton decay during preheating is bosonic 
four-legs parametric resonance, as it was assumed in the earlier papers on preheating Q. In this case, parametric 
resonance is actually amplified by the trilinear interaction, as we show in Section llVl 

To reach this important conclusion, we shall understand how preheating works with a three-legs interaction. For 
this we shall consider an idealized situation when other interactions are switched off. 

To get an insight into how preheating occurs through a trilinear vertex i<7(/>x 2 , consider the effective frequency of 
X particles caused by the inflaton field <p(t) = $ sin mi coherently oscillating around the origin 

fc 2 

w xW = -^+ml + a(j>{t) . (I) 

The bare mass of the x W1 U typically satisfy to 2 << a $. The trilinear coupling then leads to a tachyonic mass for x 
whenever tf> < 0, which happens during half of each inflaton oscillation. Correspondingly, the modes with k 2 < a \ <p\ a 2 
will be exponentially amplified during a portion of each half-period of the inflaton oscillations. Therefore preheating 
arising from a trilinear vertex shares features of both parametric resonance and tachyonic preheating, so we call 
this process tachyonic resonance. Such tachyonic resonance leads to very efficient production of x particles, so that 
preheating in this case concludes within the first few oscillations of <\> ■ 

Aspects of parametric resonance in the presence of a trilinear interaction appeared already in [f| . Tachyonic 
resonance effects may be recognised with our present understanding in a numerical plot of in the context of 
preheating after hybrid inflation. More generally, combined tachyonic and resonant effects appeared in the context of 
dynamical symmetry breaking @, and after new inflation hd]. The tachyonic effect from a trilinear interaction was 
studied in the context of the cosmological moduli problem 14] . In this paper, we systematically construct the theory 
of tachyonic resonance from trilinear interactions in preheating after chaotic inflation. 

The origin of tachyonic resonance is similar to that of parametric resonance. Without expansion of the universe 
(a(t) = 1), they both can be understood using the stability/instability chart of the Mathieu equation. Introduce 
dimensionless parameters entering the frequency (JTJ, A = q = 2% (as customary for the theory of preheating). 
The (A, q) half-plane (A > 0) is divided into stability and instability stripes (see, e.g. 0|). The line A — 2q divides 
the stability/instability chart into two regions. The case of parametric resonance for the four-legs interaction g 2 4> 2 x 2 
can be mapped onto the theory of stability/instability chart above the line A = 2q studied in The region below 
that line corresponds to what we call tachyonic resonance. On the formal side, we will extend the theory Q of 
parametric resonance "below the line A = 2q" , which requires some new technical elements. The analytic theory of 
trilinear preheating will be constructed in Section IlIII 

In this paper, we study the effects of the trilinear <f> x 2 interaction both during preheating and during the rescat- 
tering/early thermalization stage. Preheating produces a spectrum of particles predominantly in the infrared (low 
momentum) region. While 2 «-> 2 interactions can lead to kinetic equilibrium, thermal equilibrium requires an increase 
of the average energy per particle, which in turn requires particle fusion. The beginning of particle fusion is already 
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visible during rescattering, as the simulations of [9| indicate. However, this process is expected to complete on time 
scales which are much longer than the ones we can run in a lattice simulation. Thermalization after preheating has 
been studied in detail only in limited intervals of time and momentum [jj flfl ITfil | . An analytic treatment of the 
dynamics of thermalization is possible only in the simplest cases, for instance, in the pure A0 4 model |l6|. where the 
expansion of the universe can be scaled out and no mass scale appears. In an expanding universe in the more general 
case of a massive inflaton it is much harder to do analytically [flLllCj. If only quartic interactions are relevant, particle 
fusion proceeds most effectively by combining two vertices in a 4 — > 2 process. If trilinear vertices are also present with 
comparable coupling strength, 3 — > 2 processes become possible. Hence, trilinear vertices can be expected to play an 
important role in setting the thermalization time scale. Yukawa couplings can also be important for thermalisation, 
see e.g. 01 1 - I n Section [V] we present the results of our numerical simulations and we discuss the effect of the 
trilinear interaction on the dynamics after preheating. We focus in particular on the evolution of the equation of 
state. 

As we already mentioned, during preheating four-legs interactions will often dominate over three-legs interactions. 
In many cases we also may expect the presence of higher-order interactions like -p-x 2 (/) 3 , jj2X 2< t> 4 > etc. These non- 
renormalizable terms arise for example in supergravity theories. The cut-off scale M of the effective theory is expected 
to be close to the Planck scale. Planck-suppressed interactions are often supposed to lead only to perturbative 
contributions to reheating. The situation, however, may be very different for chaotic inflation because of the large 
inflaton amplitude right after inflation, $ ^ 0.1 Mp. (Clearly, one has also to make sure that non-renormalizable 
terms do not affect the inflaton potential during inflation). In fact, surprisingly, non-renormalizable terms tend to 
dominate over renormalizable terms during preheating, as long as g 2 < -A, where g 2 is the coupling of the four- legs 
interaction. 

For M close to the Planck mass, <&/M < 1 during preheating, and the effective mass squared of x is then actually 
dominated by the (fix 2 /M interaction. This vertex leads to tachyonic growth of the quanta of x, similar to what 
occurs with a three legs interaction. Note that the efficiency of this process is governed by the dimensionless parameter 
q ~ $ 3 /(m 2 M), which is very large at the end of inflation. This also opens up the possibility of very efficient preheating 
even if the inflaton does not couple to the matter sector through renormalizable interactions 2 . We consider tachyonic 
resonance due to non-renormalizable interactions in Section IVII 

We briefly summarize our results in section lYm 

II. PROPERTIES OF POTENTIALS WITH TRILINEAR INTERACTIONS 

In this section, we describe the models that we will consider in the following, and which are characterized by a 
trilinear interaction <fix 2 between the massive inflaton (f> and a light matter scalar field x ■ I n the presence of a 4>x 2 
vertex (independently of the presence of a <ft 2 x 2 interaction), we have to include a x 4 self interaction as well for the 
potential to be bounded from below, and therefore stable. This means that there are several effects of the fields 
dynamics working simultaneously. 

We consider two examples of potentials with three-legs interactions, which stress different possibilities of preheating. 

1) A simple potential which encodes the trilinear interaction between the massive inflaton <f> and another scalar 
field x is supplemented with a quartic self-interaction of x but without a four-legs interaction 4> 2 x 2 

m 2 a A , 
V= —^+2 <l3X+ 4 X (2) 

where a has the dimension of mass. Here we assume that the bare mass of x is negligible with respect to the inflaton 
mass. It is instructive to rewrite @ in the form 

^H+^ 2 ) 2+ i( A -^)* 4 - (3) 

Values A < a 2 /2m 2 are not allowed, since in this case the potential is not bounded from below. The limiting case 
A = a 2 /2m 2 is characterized by the exact flat direction <\> — —ax 2 /2m 2 where V = 0. The shape of the potential 
in this case is shown in Fig^ Finally, for A > a 2 /2m 2 the flat direction is lifted and the potential admits a single 
minimum at <f> = x = where V = 0. For numerical estimates, we will often take A = a 2 /m 2 in the following. 



1 Trilinear interactions may also play an important role for electroweak baryogenesis and leptogenesis, see llSl and references therein 

2 If a Bose condensate associated with a flat direction forms during inflation, non-renormalizable interactions may also play an important 
role during preheating, as noticed in lUi . 



FIG. 1: The potential © for A = <r 2 /2rn 2 with an exact flat direction. For A > a 2 /2m 2 the flat direction is lifted and the 
single minimum occurs at <j> = x — 0. Inflation occurs as cf> rolls down the single valley at the top (positive cj>) and preheating 
occurs as it oscillates between the single- valley and the double- valley. 



The inflationary stage is driven by the 0- field in J3J) if the effective mass of the \ ne ld is large, i.e. if a tfi(t) > H 2 > m 2 
during inflation, and the x field then stays at the minimum of its effective potential, at x = 0. After inflation (j) begins 
to oscillate around the minimum. We shall consider creation of particles of the test quantum field x due to its trilinear 
interaction with the background oscillations of the classical field <f>. 

When the inflaton passes the minimum of the potential, a 4>{t) becomes negative and x acquires a negative mass 
squared around the local maximum in the x-direction at x = (see Figure^). As a result inhomogeneous quantum 
fluctuations of x with momenta k 2 < |m^ eff (£)| = |cr(j!>(i)| grow exponentially with time, in a way similar to the 
tachyonic instability in hybrid inflation discussed in [(J. The difference here, however, is that tachyonic instability 
occurs only during the periods when the background field (f> is negative and ceases during the parts of the <f) oscillations 
when (p > 0. 

2) We will also consider the potential 

\r 171,2 jfi i a j. 2 . 9 2 ,2 2 , ^ 4 / A \ 

^ = Y? + 2« + X (4) 

with the additional four legs interaction. The potential is bounded from below for A > 0. For < A < a 2 /2m 2 , the 
potential admits two minima at non vanishing (j> and x where V < 0. We mostly focus on the case A > a 2 /2m 2 , where 
there is a single minimum at (f> = x = with V = 0. Close to this minimum the shape of the potential is similar to 
@- 

The potential (@J) may arise for instance from the superpotential 

which gives A = g 2 /2 and a = gm (for the real part of x). Preheating in the theory |0J occurs through a combination 
of effects due to the <fi 2 x 2 interaction and the (f>x 2 interaction. In the Section Hvl we discuss the range of parameters 
for which the different interactions dominate. 



III. ANALYTICAL THEORY OF TACHYONIC RESONANCE 



In this section, we study the production of quanta of a scalar field x whose effective mass squared periodically 
becomes negative in time due to its interaction with the classical inflaton field oscillating around its minimum. For 
definiteness, we consider the potential J2J- The backreaction of the x 4 self- interaction is negligible during the early 
stage of preheating, so that only the 4>x 2 term is important. 
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The quantum field x is decomposed into creation and annihilation operators with the eigenmodes Xk (t) & % ^ x , where 
k is the co-moving momentum. After the usual rescaling Xk — * a 3 ^ 2 Xk, where a(t) is the scale factor of the (spatially 
flat) FRW universe, the temporal part of the x-modes obeys the equation 



Xk+ulxk = 0, (6) 



where 



^ = ^ + <7<K*) + A, (7) 

A = —3d 2 /4a 2 — 3d/2a, and a dot denotes a derivative with respect to the proper time t. The choice of the positive 
frequency asymptotic solution in the past fixes the initial conditions while the condition XkXk ~ XkXk — i fixes the 
normalization for the solution of equation 0. In this section we consider solutions of this equation for harmonically 
oscillating (f)(t). 

Let us first neglect the expansion of the universe, a(t) = 1. The background solution for the inflaton is then given 
by 4>(t) = <& sin(mt) with a constant amplitude <f>. In this case, Eq.@ reduces to the canonical Mathieu equation 

Xk + (Ak-2qcos2z) X k = (8) 

where mt = 2z— %, A). = ^j- and q = and a prime denotes the derivative with respect to z. Right after inflation 
we have $o ~ 0.1 Mp, while m ~ 10 13 GeV in order to match with the observed CMB anisotropies. The g-parameter 
is then very large, for instance qo ~ VA10 5 >> 1 for A ~ a 2 /m 2 . 

The stability/intsability chart (A,q) of the Mathieu equation is divided into bands. We are interested in unstable 
solutions, which correspond to the amplification of the vacuum fluctuations Xfcj m the regime of large q. Let us inspect 
the effective frequency w 2 . Modes with momenta A k > 2q have positive lo\. We can view Eq.© as the Schrodinger 
equation with the periodic potential, where for A k > 2q the waves propagate above the potential barrier and are 
periodically scattered by its peaks. They are amplified at the instances when ujf, is minimal, i.e. in the vicinity of 
z = Itt which corresponds to <fi — — $ in our case, and remain adiabatic outside of those instances. The whole process 
can be considered as a series of scatterings in a parabolic potential, which approximates u> 2 around its zeros. The net 
effect corresponds to broad parametric resonance, as described in Q- 

However, for < Ak < 2q the frequency squared of the modes with k 2 < a $ becomes negative during some 
intervals of time within each period of the background oscillations. These are the modes we consider in the rest of 
this section because they give the dominant contribution to the number of produced particles. For those, we still can 
use the method of successive scatterings, but each individual scattering cannot be approximated as scattering from a 
parabolic potential. This is where we have to modify the theory. 

Suppose the inflaton oscillations begin at some initial time t = to where w 2 (io) > (take 4>(to) = Consider the 
period during the (j + l) th oscillation of the inflaton, from t = tj to t = tj+i where tj = to + 2ir j/m. The frequency 
squared of the modes with k 2 < a $ is negative during a fc-dependent interval 

n 2 (t) = ~io 2 k {t) > for t- j <t<t+ (9) 

between the "turning points" t^ - and ijy where u> 2 = 0. Now we can view Eq.10 with < Ak < 2q as the Schrodinger 
equation with a periodic potential, where the wave periodically propagates above and below the potential barrier. 

Except in the vicinity of the turning points, we have \u>k\ << \u>k\ 2 and \&k\ << \^k\ 3 (including imaginary values 
for u>k) for q >> 1, so we may use the WKB approximation to solve for JBJ. For t < ijT. (above the barrier), we have 
a superposition of positive- and negative-frequency waves 

where the integral is taken over the intervals between to and t where w 2 > 0, and oP k and f3 3 k are constant coefficients 
(in the adiabatic approximation) during this time interval, determined as we describe below. These correspond to 
the Bogoliubov coefficients (normalized as \a k \ 2 — l/3fc| 2 = 1) an d the initial vacuum for t — > to is defined by the 
positive-frequency mode, a k = 1, f3 k = 0. The adiabatic invariant 



= \Pl\ 2 (11) 



corresponds to the occupation number of the ^-particles after j inflaton oscillations and will be the major object of 
interest. 
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For t k - < t < fjy (below the barrier), the WKB approximation gives a superposition of exponentially increasing 
and decreasing solutions 

» (,) " vmw) exp (" j£ m ' )dt ') + A exp (£ nii, ' )d '') (12) 

where and b 3 k are constant coefficients (normalized as a k b k * — a k * b k — i). Finally, for t > i^, we have Xk(t) — 
Xfe +1 (^) as given by i|l(J[l in the WKB regime with the shift j —* j + 1 

The solution (|12|l with non-vanishing b k corresponds to an exponentially fast prodution of particles. The WKB 
approximation in this tachyonic (below the barrier) stage is accurate for A k < 2q — 2^fq. The duration of the tachyonic 
stage decreases with increasing k, that is with growing A k . For 2q — 2^fq < Ak < 2q there is a short tachyonic stage, 
during which however the WKB approximation does not apply. As Ak increases towards 2q, the x-modes spend less 
and less time in the tachyonic regime, which becomes less and less efficient. For Ak > 2q — 2Jq, particle production 
occurs in a small interval around cf> = — providing a smooth limit with the case of broad parametric resonance for 
A k > 2q. 

To evaluate the growth of rik during one oscillation of the inflaton, we have to match between the successive 
approximate solutions l|lUfl and l|12|) around the turning points, where the WKB approximation is inapplicable. This 
may be done by formally extending the domain of definition of \k to the complex plane, see Appendix A. Except 
for the normalization conditions, the procedure is similar to the calculation, in quantum mechanics, of the (spatial) 
evolution of a wave function in the quasi-classical regime |2iJ. We then find the "transfer matrix" between the 
Bogoliubov coefficients after and before the (j + l) th oscillation of the inflaton 3 

^ ( 1 



ft 



■i+i 




^ .... I (13) 



where 



x{ = / n k {t')d£ (14) 



and 9 k is the total phase accumulated from to to t k j during the intervals where cu k > 0, 9 k — J^ 3 dtu>k(t) 

When the expansion of 
case 9 k = 8 k + j Qk where 



When the expansion of the universe is neglected, and Xk do not depend on j. We have furthermore in this 



e> k = / u) k (t')dt' (15) 

is the phase accumulated during one inflaton oscillation when uj k > 0. It is then easy to apply H13[) (j + 1) times 
recursively. With the initial conditions a k = 1 and f3 k — 0, one finds the occupation number of the x-particles in the 
fc-mode (with k 2 < a $) after j oscillations of the inflaton to be 

n{ = \Pi\ 2 = exp(2jX k ) (2 cosefc) 2 ^" 1 ) . (16) 

This simple formula is the main analytic result of our paper, exp (2Xk) gives the occupation number after the first 
oscillation. For the trilinear interaction (0), one finds from (|14|) and H15[) . in terms of the variables in the Mathicu 
equation ijHJ 

X k = J + L n k (z) dz = 2 ^/2q - A k E [z k ■ ^ ^ (17) 



3 This may be expressed as usual in terms of reflection and transmission coefficients satisfying |-Rfc| 2 + |£*fc| 2 = 1, up to exponentially 
small terms, oc exp I — X^j , which are beyond the accuracy of the WKB approximation. 
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and 



(18) 



where £?(# ; m) denotes the incomplete elliptic integral of the second kind with amplitude 9 < ^ and parameter m 
1151 . Here we have defined 



~z h =\ arccos (^j e [0, J 

so that the turning points are given by = 7r ± z k ■ 

The functional dependence of X k , which controls the efficiency of particle production, is 

X k = 2^2q~f(^- 



(19) 



(20) 



In the interval we are interested (A k < 2q), we found a good approximation / (y) ~ 0.6 (1 — y) (verified numerically). 
Therefore, we can use the following accurate approximation 



A k + 2xy/q 



(21) 



with T - r (3/4) 
WlWl X - 2V2 r(5/4) 



0.85. 



The analytic formula l|16[) derived with the WKB method, gives a pretty good approximation to the actual field 
dynamics. To show this, we plot in Figure the occupation number n° k in 1)16(1 as a function of A k , after j ' = 1 and 
j = 4 oscillations, for the value of the parameter q = 20. Numerical calculations (dots) coincide with the results of 
the analytical curve. 



Log [n_kl ] 




Log [n_k4 
60 



A k 




«!/• A k 



FIG. 2: logn^. derived from 1161 as a function of A k , after j = 1 (left) and j = 4 (right) oscillations of the inflaton, for q — 20. 
The dots correspond to numerical solutions. 



Quite interestingly, after several inflaton oscillations (j > 1), there are modes with negative frequency squared 
for which the tachyonic growth does not occur (see right panel of Figure El • These form stability bands, which 
shrink to lines for q >> 1, given by cosOfc = (see 11161 ). In some sense, the situation is reversed with respect 
to the conventional parametric resonance for q -C 1, where only some values of k, forming the instability bands, 
"resonate" with the frequency of the inflaton and undergo an exponential growth. Here, all the modes are amplified 
by the tachyonic instability except the ones belonging to the very narrow stability bands. These are described in more 
details in Appendix B. Beyond the WKB approximation, each stability line is split after each inflaton oscillation. After 
several oscillations, the net effect is the formation of stability bands of finite width 5A k ~ y/qe~ x ^. As q decreases, 
these bands become wider, while the distance between them decreases as AA k ~ y/q (see appendix B) ; at the same 
time the efficiency of particle production decreases as X k ~ Jq from When q < 1, the width of the stability 

and instability bands are comparable, and the tachyonic effect becomes indistinguishable from (narrow) parametric 
resonance (as is clear from the stability/instability chart of the Mathieu equation in this region of parameters). 
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The expansion of the universe can be incorporated "adiabatically" , which is to say that the above formulas will stay 
the same but the parameters will depend on a(t). In an expanding universe, we may essentially repeat the discussion 
leading to without any modification, because of the hierarchy of scales $/$ ~ vA ~ H « m « V <7<I> at 
the beginning of preheating. The squares of the physical momenta now redshift as a~ 2 with the scale factor, which 
is faster than the amplitude of the inflaton oscillations, $ ~ a -3 / 2 . This means that more and more modes satisfy 
k 2 /a 2 < er$(i) and are therefore amplified by the tachyonic instability. The modes which are amplified latest give 
the highest contribution to the comoving phase space volume. However, this effect is subdominant with respect to 
the decrease with time of the rate of particle production X^. Using (|21|) . X 3 k during the j th inflaton oscillation may 
be estimated as 



-, i n ' 



go a a } 

where go is the initial g-parameter and aj is the mean value of the scale factor during the j th period of particle 

production. Since the second term dominates this expression, X 3 k decreases approximately as Xj, ~ ^fq ~ 3 ^ 4 with 

the expansion. As in the case of broad parametric resonance, the phase 6 k in (| 1 31) is random in an expanding universe. 
This results here in broader "stability bands" as compared to flat spacetime. Contrary to preheating due to the cj) 2 x 2 
interaction, the interference term in 9 k in (|13|) is never large enough to decrease the total number of particles during 
some oscillations of the inflaton. In conclusion, the main effect of the expansion of the universe is a decrease of the 
efficiency of particle production at each inflaton oscillation, due to the decrease of the effective g-parameter. This 
tends to lengthen the preheating stage. If the g-parameter decreases up to q ~ 1, we are then in the regime of narrow 
parametric resonance, which is shut off by the expansion of the universe as discussed in In general, however, 
the actual dynamics at the end of preheating is rather complicate and includes backreaction of the quanta produced, 
rescattering and self-interaction. 

The analytic derivation leading to the formula Qlfijl is valid for tachyonic values of uj 2 , which can include 4>'\ 2 
and 4> 2 x 2 interactions, as far as the three-legs coupling dominates. The tachyonic effect is stronger for stronger 
three-legs couplings (relatively to the four-legs couplings). In the limiting case, when the 4> 2 x 2 interaction is absent, 
preheating is so efficient that it completes typically after the first or second oscillation of the inflaton homogeneous 
mode. Indeed, as the x-modes are amplified, their variance (x 2 ) grows, and the contribution of the A x 4 self-interaction 
to their frequency squared, 3A(x 2 ), eventually becomes significant. The tachyonic instability is then shut off when 
(X 2 ) ~ ct$(*)/(3A). (In principle, some modes may still be amplified by parametric resonance during the next inflaton 
oscillation.) Let us roughly estimate the value of the g-parameter for which this occurs after the first oscillation of 
the inflaton. For this me may neglect the expansion of the universe and work with the value of q during the first 
"tachyonic interval" , which is about one half of its initial value go at the end of inflation. From H21(l , the variance of 
X after the first oscillation of the inflaton may be estimated as 

For instance for the value A = cr 2 /m 2 , this is equal to tr$/(3A) for g ~ 50 and go = 100. The corresponding energy 
density in x-particles may be estimated in the same way 

which is about 5 percent of the initial energy density of the system, pi ~ i m 2 $ 2 . These estimations are over- 
simplified but the orders of magnitude are correct. In particular, most of the energy density is still in the coherent 
inflaton mode at the end of preheating, and it will mainly decay during the "rescattering" stage. 

In this section we outlined in detail the formalism of particle creation from a <7^>x 2 interaction. The formalism in 
fact is more general, and can be applied to the more general case when both ct0x 2 an( i 5 2 ^ 2 X 2 interactions contribute. 
Instead of repeating the procedure for this case, in the next section we will discuss when each of them may dominate. 



IV. COMPARISON OF THE EFFECTS OF THREE AND FOUR LEGS INTERACTIONS 



In this section we will address the important question: Which interaction, three- or four-legs, is dominant during 
preheating ? We shall estimate the parameter ranges for which parametric resonance or tachyonic amplification will 
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dominate at the stage of preheating, and discuss how the three-legs interaction enhances parametric resonance even 
if the trilinear term is subdominant. 

After inflation the inflaton oscillates sinusoidally with a decaying amplitude <&(t) = <I>o/a 3 / 2 . When the backreaction 
of the x 4 self-interaction in Q is still negligible, the time-dependant frequency of the x-modes is given by 

tol(t) = k 2 + cr sin(rat) + g 2 <P(t) 2 sin 2 (mi) , (25) 

where k denotes the physical momenta. Without expansion of the universe, this reduces to a particular case of the 
Hill equation 

Xk + [ A k - 23 cos z + y cos 2 Xk = , (26) 

where its parameters are slightly different from those in ©, namely z = mt + tt/2, = k 2 /m + 94/2 and we have 
defined the dimensionless quantities 

cr $ o 2 $ 2 
c/ 3 = — r and c/ 4 = 5- ■ (27) 

771 777 z 

Equation (|26[) is an equation with periodic coefficients, therefore, according to the Floquet theorem, it admits 
exponentially unstable solutions. However, the structure of the stability/instability chart will be denser and more 
complicated than that of the Mathieu equation. In any case the expansion of the universe for the interesting case 
of large parameters c/3 and c/4 makes the stability/instability chart relevant only heuristically. Therefore we restrict 
ourself to qualitative arguments based on the physical picture of particle creation from three- and four- legs vertices. 

First consider the case c/3 < c/4, for which the frequency squared (|25|l is negative in the interval — C/3/C/4 < sin(mt) < 0. 

J 2 a 2 

This leads to a tachyonic amplification of modes with < j^. In addition to this effect, modes are produced by the 
four-legs interaction when the non adiabaticity condition, \tOk\ > |wfe| 2 j is satisfied. This happens in small intervals 
Asin(m£) ~ 2/q^ 4 around sin(mi) = — C/3/2C/4, as may be found by applying the calculation of 3]. The typical 
momentum of the modes amplified by this process is given by ^ = + Inspecting the right hand side of this 

expression, we can break the range of c/3 and c/4 into two regions, c/3 < q^ 4 and c/3 > c/4^ 4 . 
For the case 

93 < qT (28) 

the periods of tachyonic instability are negligible, and particle production occurs due to the non-adiabaticity in very 
small intervals mAt << 1 due to the four- legs vertex. Yet, the presence of a three- legs vertex facilitates the effect. 
Indeed, in this case, we may use the method of successive parabolic scatterings to evaluate the growth of the 
occupation number during one event of particle production. After tedious calculations we find 



2 e"™ 2 ( 1 - sin 9 \/e™ 2 + 1 ) 77^ (29) 



77 i+ l = 1, 

where 9 is the phase of waves accumulated between successive scatterings. In an expanding universe and for large 
q parameters, 9 can be treated as the random phase so that the resonant effect is stochastic. (In a sense the 
stability/instability bands are smeared, see [3] for details.) In formula l|29(l 

4c/4 /2 

which is different from k 2 = q 4 for the case when only the four legs interaction is present. Thus, for the range 
of parameters l|28|l when preheating is dominated by parametric resonance of the four-legs interaction, the trilinear 
interaction enhances particle production due to the additional factor e' Rq3 ^ 4q4 in l|29|l . 
In the opposite case 

93 > qT (31) 

the tachyonic effect is dominant, and we have to use formula (|I6|I of the previous section. 

We are especially interested in the SUSY case. Let us take the superpotential which leads to c/3 = Jqi- For 
c/3 , c/4 ^ 1 this corresponds to the case Q28H ■ This means that in the SUSY theory where both three and four legs 
interactions are present, the stage of preheating is dominated by the four legs interaction. 
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Suppose that initially we have condition (|28|l . as in the SUSY theory. Consider q^ 4 jqz as function of time. The 
inflaton amplitude redshifts as $ ~ (mt)~ 4 with the expansion of the universe, so that the critical ratio q Z J 4 jq± 
decreases with time as (mt)^ 1 ! 2 . If initially this ratio is larger than unity, it decreases with time and later on the 

situation reverses, g 3 / 4 /q± becomes less than unity, and the three legs vertex dominates. 

We conclude that in many cases, notably in the SUSY theory, the stage of preheating is dominated by the effects of 
the four legs interaction, while the later stages are dominated by the three legs interaction. As we will see in the next 
section, independently of the dominant mechanism for preheating, the trilinear interaction dominates the following 
turbulent dynamics. 

V. NON-LINEAR DYNAMICS AFTER PREHEATING WITH A THREE LEGS INTERACTION 

So far we have analytically studied preheating with a three-legs interaction and compared the effects of three- and 
four-legs interactions during the linear stage. After the system becomes non-linear, the analytic theory is replaced by 
numerical simulations of the dynamics. In this section we report the results of those simulations. 

We will consider three different models, one with only a four-legs interaction, one with only a three-legs interaction, 
and one with both three- and four-legs interactions. Let us define these three models as 



Model A: 


qu 


= 10 4 


q X i 


= 5 x 10 3 


<73i 


= 


(32) 


Model B: 


q±i 


= 10 4 


qxi 


= 5 x 10 3 




= 10 2 


(33) 


Model C: 




= 


j qxi 


= 10 4 


q3i 


= 10 2 


(34) 



(35) 

where we have introduced the dimensionless quantities q^ = g 2 ^\jm 2 , q 3i = <j& /m 2 and <f> — 0.193 Mp is the 
initial amplitude of the inflaton zero- mode at the end of inflation. The parameter q xi = AC^/m 2 characterizes the 
self-interaction of the x-field. 

Model A is a standard preheating model with parametric resonance due to a four-legs interaction, without three- 
legs coupling. Model B contains both three- and four-legs couplings and corresponds to the SUSY case discussed in 
Section^ with X = g 2 /2 and cr = gm in @. Model C corresponds to the potential (J2J without cf> 2 % 2 interaction (when 
the exact flat direction shown in Figure^is slightly lifted, A > ct 2 /(2to 2 )) and thus allows us to see the pure effects 
of tachyonic resonance due to the three-legs coupling. A \ 4 self-interaction is present in all models (to ensure the 
stability of the potential in the presence of the cj> x 2 interaction, see Section [nj . We will see that this self-interaction 
does not change qualitatively the dynamics. 

We used the LATTICEEASY program [2l| to perform numerical lattice simulations for the theory Q in an 
expanding universe. For the three models, we used a 3-dimensional lattice with N — 256 3 points and a comoving edge 
size L — 10/m. This allows us to probe comoving momenta in the range 0.6 m < k < 140 m. Here we are interested in 
the behavior of the system both during and after preheating. Long time integrations require both a high ultra-violet 
(UV) cutoff for the modes and a high resolution of modes in the infra-red (IR) part of the spectra. This motivates 
the value of q<u chosen, which is large enough to produce efficient preheating, but also small enough to produce quick 
rescattering [l(| ■ For the three models, the simulations were performed up to a total period of m At = 1000 starting 
from the end of inflation. We also checked the robustness of our results by performing 2-dimensional simulations with 
both a larger UV cutoff (smaller lattice spacing) and a higher resolution in the IR (larger box size). 

As output of the calculations, we present the evolution of the mean values of the fields, the evolution of the 
occupation number spectra, and most importantly the evolution of the equation of state. For the definitions of these 
quantities, see e.g. [lOf. As we will see, after the preheating stage model B (with both types of interactions) will 
behave qualitatively similarly to model C (with only a three legs interaction), but significantly differently from model 
A (with only a four legs interaction). We can thus conclude that the three legs interaction dominates the dynamics 
of model B during this stage. 

The evolution with time of the homogeneous inflaton zero- mode in the 3 models is shown in Figure |3] In model 
A without the trilinear term the inflaton mean decays to zero, but in the presence of the trilinear term in models B 
and C the backreaction of the quanta produced displaces the global minimum of the potential away from <j) = \ = 0- 
Therefore (f> develops a non vanishing VEV (j) — 4> ~ S(f>, which is in very good agreement with the one found from the 
Hartree approximation <f> — — § (x 2 ) / ( m2 + 9 2 (x 2 )) where (% 2 ) is the variance of \- 

The occupation numbers of x an d S<t> quanta are shown in Figures El and El for the three models. We can 
distinguish three qualitatively different stages: linear preheating, the short violent non- linear stage and the long 
turbulent stage where spectra cascade towards saturated distributions. The preheating stage in which the spectrum 
of x grows rapidly occurs faster in model B than in model A, showing the effect of the trilinear term on parametric 
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FIG. 3: Mean value of <f>, in units of $o and multiplied by a 3 ^ 2 (t), as a function of mt, for models A (left), B (middle), and C 
(right). 



resonance. The theory of this effect (impact of three-legs vertices on parametric resonance) was given in Section IV. 
Meanwhile in model C, where preheating occurs through tachyonic resonance only, this stage takes only one inflaton 
oscillation, showing the efficiency of the effect. 



Number: x Number: x Number: x 
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FIG. 4: Occupation number n£ as a function of k/m for models A (left), B (middle), and C (right). The spectra go from early 
times (red plots) to late times (blue plots) with spacing At = 10/m. 
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FIG. 5: The same as in Figure |1] but for S<f> quanta. 

After preheating ends there is a brief violent stage of rescattering, followed by a much longer stage of weak turbulence 
during which the spectra cascade towards UV and IR modes, gradually approaching kinetic equilibrium. We can see 
in Figures 01 and that this cascade occurs faster in the presence of a trilinear term, i.e. in models B and C compared 
to model A. In Figure we plot k^io^n^ for 5<fi quanta, which corresponds to the contribution of these modes to the 
total energy density of the system at the end of the simulations. We see that in the presence of the trilinear interaction 
the spectrum is spread much more strongly towards the UV. This means that more 8(f) quanta are relativistic in this 
case. This will be important for the following discussion. 

Now we turn to the evolution of the equation of state (EOS), see Figure [7] In model A, where the trilinear 
interaction is absent, the equation of state rises to roughly 1/4 and then starts falling back towards zero. In this 
model the EOS will never reach w = 1/3 because the massive inflaton particles cannot decay completely 01 ■ I n the 
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FIG. 6: Spectra k 2 tut n,k/m 3 as a function of k/m for S(f> quanta at rat = 1000. The spectrum peaked in the IR is for model 
A while the spectra for model B (above) and model C (lower) are spread out much further in the UV. 



presence of a trilinear term, by contrast, the equation of state jumps to a plateau value and does not decrease. In 
our simulations of classical fields this value is slightly below 1/3. This is the most important reason to introduce a 
three-legs interactions. The trilinear term dramatically affects the evolution of the equation of state. It is expected 
that much later in the evolution, when quantum effects become important, the decay of the massive inflaton will 
result in an asymptotic radiation EOS. 
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FIG. 7: Equation of state w = ptot/ptot of the system, as a function of mt, for models A (left), B (middle), and C (right). 



To understand the behavior of the EOS, we investigate how the fraction of relativistic particles is evolving in the 
system. Let us consider the fraction of quanta whose physical momenta k/a are larger than their effective masses 



Ik 



^■rel Jfc>am c ff 

n 



d 3 k 



Ilk 



J d 3 knfe 



The effective masses are given by 



m |,eff 



s (t) = 5 V + ^ + .9W}+3A(x^> 



m + g A (x 2 ) 

-2 12 



(36) 



(37) 
(38) 



for 4> and x respectively. We calculate and plot in Figure 03 the fraction of relativistic particles. 

After the preheating stage, the effective mass of the %-quanta is dominated by the term 3A (x 2 ), which is diluted with 
the expansion of the universe as, approximatively, m XjC ff ~ yO?) ~ This is the same rate as the redshift 

of physical momenta, so that the number of relativistic x~°L u anta increases due to the flow of comoving momenta 
towards UV modes. On the other hand, the effective mass of the 50-quanta is dominated by the constant bare mass 
to. Therefore, these quanta become less relativistic as their physical momenta are redshifted with the expansion, 
except if the flow of comoving momenta towards UV modes is as fast as a(t). We see from Figure |H1 that this only 
occurs in the presence of a trilinear term. In short, model A moves towards matter domination because the cf> particles 
cannot decay completely and become increasingly non-relativistic. In the presence of a trilinear interaction, however, 
the <f> particles become increasingly relativistic. The inflaton field should decay completely on a much longer time 
scale through the quantum process 4> — > XX- When the occupation numbers have been diluted by the expansion of 
the universe, n/j/a 3 << 1, this process should dominate the dynamics of model B and C. In this regime, the typical 
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FIG. 8: The fraction of quanta that are relativistic, 2^ as a function of mt for Sep (dark) and x (g re y) i n models A (left), B 
(middle), and C (right). 



time of infiaton particle decay is given by t ~ 1/T ~ 87rm/er 2 ~ 10 8 /m for — 100. This stage of the evolution 
cannot be captured by simulations of the classical fields dynamics, however. 



VI. PREHEATING WITH NON-RENORMALIZABLE INTERACTIONS 

In this section we study preheating due to the presence of non-renormalizable terms in the theory 

fx 2 



AC = -A„ 



M 



n-1 



(39) 



These terms may occur alongside the usual three- and four-legs vertices, or in cases where the only interactions are 
the non-renomalizable ones. We show that, quite surprisingly, the non-renormalizable terms 1|39[) alone lead to very 
efficient preheating. We also show that, even in the presence of renormalizable interactions, non-renormalizable ones 
tend to dominate the early stage of preheating. 

If M is close to the Planck mass, the amplitude of the infiaton satisfies $/Af < 1 during preheating and the theory 
is in a controllable regime. If all the couplings A„ are of order one, the dominant interaction in Ij39(l is then the one 
with n — 3. Here we will consider the more general case where n is odd. For n even, there is no tachyonic amplification 
and preheating occurs due to parametric resonance. 

Let us begin with the theory of x particle creation from the oscillating infiaton due to the interaction (|39|) with n 
odd, without expansion of the universe. The Xk eigenmodes obey the equation 



where mt = 2z — -| , Ak 



n ~ rn 2 M 



X'i + (A k - 2q n cos" 2z) X k = 0, (40) 

2A„*" p or M _ Mpj Xn _ and m _ 1Q 13 Q eVj we have Qn _ 10 12/10« 



and q 

at the beginning of preheating 

The method of solving this equation for q n ^S> 1 will be based on a generalization of the method of Section ITTTI which 
we developed for equation JSJl for the regime of tachyonic resonance. 

Repeating the analysis of section ITTTI and the Appendices, we find that the behavior of occupation numbers of % 
particles with time is given by the formula which generalizes (|16ll 



nf\k) = exp(2jXl n) ) (2 co S e£ 



(n) 



20-1) 



with 



and 



X 



(n) 



\j2q n cos™ 2z — Ak dz 



(41) 



(42) 



e 



V 'A k - 2q n cos™ 2z dz, (43) 
where Zk tU — ^arccos ( J ■ The condition cosOjj™' 1 = defines the median locations of the stability bands. 
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To estimate the rate of particle production, we shall estimate the function -Xl given by the integral (|42|l . It turns 
out to be well approximated by the formula 



^r(^) / (a 



2 r(l + f) I \2q n 



(44) 



This reduces to 1)21)1 for particular case of n = 1. 

For illustration, we apply the obtained results for a particular non-renormalizable interaction A3 9 ^ , corresponding 
to n = 3 in formulae above. The spectrum of x particles created from tachyonic resonance in this case is shown in 
Figure El for q% = 20 and q$ = 200, after j = 4 inflaton oscillations. 



Log[n_k4] Log[n_k4] 




FIG. 9: Occupation number log n A (k) as a function of A%, resulting from tachyonic resonance due to the interaction ^ * 
after j = 4 inflaton oscillations, for qs — 20 (left) and q$ — 200 (right). The curve is obtained from the analytical formula 1411 
while the dots correspond to numerical solutions of 1401 with positive frequency initial conditions. 



Comparing to Figure |2 for n = 1, we see that, for the same not very large value of the g-parameter (q = 20, left 
panel of Figure (SJ, the WKB approximation is less accurate for n = 3 (and more generally for n > 1) than for n = 1. 
This approximation improves with the growth of q. Recall that for n = 1 the approximation works perfectly even for 
moderate q. 

To understand the subtle difference between the n = 1 and n > 1 cases, we shall estimate the range of validity 
of the WKB approximation for arbitrary n. Indeed, applying (|A2|) to 1)41) Jl. we derive the range of validity for the 
matching procedure leading to 141)) (see Appendix A) 



("-!) 
16 q n 



n/(n+2) 



Av I n 

« < 1 -\ h 45 
2q n V In 



The upper bound was already obtained in Section lllll in the particular case n = 1. However, for n 7^ 1, we now have 
a lower bound too. This explains the mismatch between the analytical and numerical results for small in Figure^ 
for 53 = 20 (left panel), which was absent for n = 1. The reason why <|41|l is less accurate for 4 ~ when n > 1 
is that, in this case, the frequency squared in (|40() varies as \z — zq\ u in the vicinity of the turning point zq, whereas 
the matching procedure used relies on the linear approximation i|Al[) . which breaks down for A^ « (see Appendix 
A for details). On the other hand, we see from Q45JI that the range of Ak for which 1)41(1 is accurate increases with q. 
This is confirmed in the right panel of Figure E| At the beginning of preheating, the g-parameter is very large and 
the analytical approximation l|41|l is very accurate. 

We see that the impact of non-renormalizable interactions on preheating can be very strong. 

Let us now include the expansion of the universe. We can repeat the previous calculations just by "adiabatic" 
scaling of parameters with the scale factor 



For non-renormalizable interactions, the g-parameter decreases with time faster for bigger n, and in particular 
faster than for the trilinear interaction. For n > 1, it decreases faster than the square of the physical momenta, 
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which means that in this case fewer and fewer modes are amplified by the tachyonic effect. The overall effect of the 
expansion of the universe is thus to decrease the effectiveness of tachyonic resonance (see Q44|)'l. and this occurs faster 
for non-renormalizable interactions than for the trilinear one. In particular, depending on the interactions which 
stabilize the potential, the preheating stage may now terminate in the regime of narrow parametric resonance, q n ~ 1 
(see Section Ijlllt ). For example, for n = 9 and q n ~ 10 12 /10 rl at the beginning of preheating, this happens after only 
2 oscillations of the inflaton field. 

Let us now consider other terms in the potential of the effective field theory describing the interactions between the 

in 2 

inflaton and the matter field X: such as a four- legs interaction and non-renormalizable terms A K j^„~2 with different 
n. The effective mass squared for x is given by 

As we already discussed, the contribution of the renormalizable interactions is generally dominated by the term 
g 2 <j> 2 X 2 m the potential. At the beginning of preheating, the amplitude of the inflaton is given by $ ~ 0.1 Mp. For 
g 2 < 4|> < 1 and all the couplings Aj of order one, we see that preheating is dominated by the term A3 9 ^ in the 
potential. The important conclusion is that it dominates over the A4 9 ^ term, the A5 9 ^ term, and the other 
non-renormalizable terms, as well as over the renormalizable terms such as g 2 (f> 2 x 2 ■ 



VII. CONCLUSIONS 



The main purpose of this paper is the study of the impact of bosonic trilinear interactions on the decay of the 
inflaton after chaotic inflation, and on the following thermalization stage. The starting motivation is that terms of 
the form (fix 11 (where <f> is the inflaton, while x represents some matter field) are necessary for a complete decay of the 
inflaton, and the trilinear interaction <f>x 2 is certainly the most common of such interactions. In particular, it appears 
naturally in supersymmetric theories. In our analysis, we found that a trilinear coupling can lead to other relevant 
effects. 

One of them is the appearance of a qualitatively new type of preheating, for which we developed an analytical 
theory. On the formal side, we extended the theory of broad parametric resonance, developed in [j| for the region 
A > 2q, to the remaining portion of parameters where < A < 2q. Preheating in this region can be best understood 
as a combination of parametric resonance and tachyonic effects, so we called it "tachyonic resonance." For half of 
its oscillation the inflaton is negative and 4>x 2 provides a tachyonic effective mass for the scalar x ■ This tachyonic 
instability recurs periodically in time, which results in the presence of resonance/stability bands, which is typical 
of parametric resonance. Depending on the relative strengths of the interactions, this effect can be dominant or 
subdominant relative to the more studied parametric resonance from a quartic 4> 2 x 2 interaction. If the trilinear 
interaction comes from a simple superpotential, the quartic interaction is dominant at early times. Nevertheless, the 
trilinear interaction tends to enhance the parametric resonance driven by the four-legs interaction, due to the fact 
that the resonance occurs when (f> « , where the relative contribution of the cubic over the quartic interaction is 
maximal. Even if initially subdominant, the cubic interaction eventually comes to dominate over the quartic one due 
to the decrease of <j) as the universe expands. Eventually, it leads to a complete decay of the inflaton, while the quartic 
interaction freezes out. 

A different noticeable effect of the trilinear term is a speed up of the thermalization of the system. The distributions 
produced at preheating are largely populated at (relatively) low momenta, and they are much more enhanced in the 
infrared with respect to a thermal distribution. Therefore, thermalization proceeds through particle fusion. As can be 
expected, the presence of a trilinear interaction enhances this effect. We confirmed this by evolving fields on the lattice 
and comparing models with and without the <f>x 2 term. If only the quartic (f) 2 x 2 interaction is present, preheating is 
followed by a turbulent rescattering stage, where quanta of <f> are also excited, and by a much longer phase where the 
distributions very slowly evolve towards the ultraviolet. In the meantime, the equation of state (EOS) of the system 
rises to an intermediate value between the one of matter and radiation, and it then evolves back towards the one of 
matter, due to the increasing importance of nonrclativistic quanta of the massive inflaton. On the contrary, in the 
cases in which the trilinear term was present we noticed a much quicker population of the UV modes, particularly 
for the inflaton field. As a consequence, the EOS remains fixed at the intermediate level throughout our simulation, 
while the fraction of relativistic quanta of 4> continues to increase. Although thermalization completes on a much 
longer timescale than we can simulate, details of this early stage can be extremely relevant for the production of heavy 
particles, gravitational relics, and for modulated perturbations, as we have described in |To| . 

Tachyonic resonance (with small variations of details) also arises for non-renormalizable interactions in <f> n x 2 with 
11 > 1 odd, and we extended our analytical formalism to such couplings. Non-renormalizable interactions naturally 
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arise in supergravity models, where they are suppressed by powers of the Planck mass. In chaotic inflation, <f> ~ 
O.IMp at the beginning of reheating, so the higher order operators are progressively more and more suppressed (this 
guarantees that the computation is under control). However, quite surprisingly, we found that a \<j) 3 x 2 /M p term 
(with A of order one) dominates over the quartic g 2 4> 2 X 2 interaction at the beginning of preheating, for the values 
of coupling g 2 which are typically considered. Such a A0 3 x 2 /M p term leads also to much stronger preheating, via 
tachyonic resonance. This opens a qualitatively new possibility for inflaton decay, even if it interacts only very weakly 
with the Standard Model sector through usual renormalizable interactions. 
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APPENDIX A: DERIVATION OF THE TRANSFER MATRIX (Q3j| FOR TACHYONIC RESONANCE 



In this appendix, we outline the procedure to derive the transfer matrix (|13[1 and we discuss its range of validity. 



First we want to follow the evolution of the modes from an oscillatory form (|10f) to an exponential form 1|12|) in the 
vicinity of the turning point t = t^- . Let us decompose the frequency to linear order 



= -sH 



dz 



(Al) 



Along the real axis of t we can only use the WKB approximation far away from the turning point. However, we 
can implement here the following trick known in quantum mechanics |20| . By analytical continuation of Xkif) to the 
complex plane, we may pass around the turning point along a path in the complex plane t which is sufficiently far 
from it so that the WKB approximation remains valid, but also sufficiently close to it so that we may use I|A1|I. It 
follows that 



dz 



1-1/3 , 2 



d 2 LO 



dz 2 



(A2) 



has to be satisfied along the contour. For the three- legs interaction |JSJ, such a contour always exists for the modes 
with Ak < 2q — which is in any case required for the approximation (|12f> to be valid (see section 3). The 

constraint (|A2|) is actually more severe in the case of non-renormalizable interactions, as discussed in Section IVII 




Ret 



FIG. 10: For sufficiently large p, the WKB approximation is applicable everywhere along both contours. 

The contour may be chosen in the lower or in the upper half of the complex plane as shown in Figure ITU1 Using the 
representation t^- — t = pe 1 ^ and varying ip from to ir, one goes from negative to positive t — t^- along a contour 
in the lower half of the complex plane. This gives 

(*«-*) 1/a - (t-t^) 1/2 e^ 2 (A3) 

so that (|10f> is matched to 

a{e- l ^ + *'^ I*- «""*(*) 
Xi - e *« (A4) 
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where we neglected an exponentially small term proportional to j3 3 k , which is beyond the accuracy of the WKB 
approximation. With our convention, 9 k is the total phase accumulated from to to £7y. during the intervals where 

i — 

io\ > 0, 9 k = f t !° J dtuik{t). Similarly, going along a contour in the upper half of the complex plane gives 

M e i(ei+K/4) p dt n k (t) 
X{ -» e *w (A5) 

Comparing (|A4|) and HA5|> to $F2$ , one finds 

K = < e-^i+^^ + ^ k eWi+"/V (A6) 



which relates the evolution of the modes before and after the turning point t — t k ^ . Repeating the same procedure 

"hi " 

(3 k +1 = e x i e -^+-/ 4 ) b{ and of 1 = e x * e ^+-/4) y (A7) 



around the turning point t — t~t- gives 



where X k is defined in (|14|l . Again, the matching procedure is accurate if the contour may be chosen such that a 
relation similar to l|A2() is satisfied along it. The transfer matrix i|13|) follows directly form (| A6|> . (|A7|) . 

A similar result can be obtained if we substitute the linear form (|A1|) in the wave equation © and observe that 
its solutions are given by Airy functions. Matching asymptotics of the Airy functions with the oscillating IjlOfl and 
exponential forms (|12fl gives the transfer matrix 

APPENDIX B: COMMENTS ON STABILITY BANDS 

In the regime of tachyonic resonance in the region where A < 2q the stability bands are very narrow for <; > 1. The 
situation is reversed with respect to conventional parametric resonance for A > 2q where for g « 1 the instability 
bands are very narrow. Formally, there is an interesting analogy between time-dependent solutions of the wave 
equations with periodic frequency and the one-dimensional Schrodinger equation for x-dependent wave function in a 
spatially periodic potential. Solving the eigenvalue problem for the Schrodinger equation with a periodic potential 
results in finding the bound states and their energy levels. Each energy level is slightly washed out due to the tunneling 
through the barriers between minima of the potential. Those energy bands exactly corresponds to the stability bands 
in the theory of time-dependent tachyonic resonance. 

Let us first consider the Schrodinger equation with the potential V(z) = 2qcos2z for < z < it, so that it has 
only one minimum. In the WKB approximation the bound energy levels are defined with the Bohr-Zommerfeld 
quantization formula Ok(^-k, <?) = 7r(n+ 1/2), where 0& is the phase of the wave function and n is an integer number. 
In our case we have 

e fc (A fc , g ) i r — _ i 

= — / yAk — 2q cos 2z dz = n-\ — (Bl) 

K n 'iarccos^ " 2 

In the time-dependent problem of tachyonic resonance this exactly corresponds to the stability band around k where 
cos6fc(z) = 0, see 

The function 8fc(Afc, q) was calculated in l|18|) . To find the position of the stability bands in the (A/., q) plane, we 
have to invert ijBlfl . We notice, by direct differentiation, that 

de k (A k: g) T 

dA k - 4 (B2) 

where T is the classical period of motion between the turning points in the potential, 

T = 2 / A a . (B3) 

7i arccos £i y<A k - 2qcos2z 

Using properties of elliptic integrals, this reduces to 
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where K is the complete elliptic integral of the first kind, which has a logarithmic divergence when its argument tends 
to 1 (see We find the accurate approximation in the range < A k < 2q 



T 



3.71 



Integrating (|B2(1 then gives 



e k (A k ,q) = 



1 - 0.235 In 1 



,A k 
o — 
2q 



Ak 

2q 



A k \ , / A k 

1 - In 1 

2q) \ 2q 



(B5) 



(B6) 



A very good approximation is obtained for a = 1.69, b = 2.31 and c = 0.46. Using the asymptotic behavior of l|B6|) 
for Ak/2q — > and A k /2q — ► 1, we find from I|B1|) the distance between the adjacent stability bands: 

AA k =d^q- (B7) 

where d ~ 3.41 for A k /2q — > and d ~ 2.73 for A k /2q — > 1. The stability bands are approximately equidistant in 
both limits. Eq. (|B7Jl agrees with Fig. [3 (right panel). 

Using the analogy with quantum mechanics, we can go beyond the WKB approximation. In particular, one can 
take the periodicity of the quantum mechanical potential into account [22| and find the exact wave function of the 
particle Xk, which will have the form of a Bloch wave function. It turns out that the only effect of the periodicity of 
the potential on the spectrum is the broadening of levels IjBljl . The characteristic width is 



SA k 



D 



(B8) 



where D — exp(— X k ) is the transmission coefficient and the function X k was calculated in l|17l I21JI . It follows that, 
to first order in A k /2q, the width of the bands varies with q as SA k ~ ^e" 1 ^, where x ~ 0.85. 

The stability/instability chart of the Mathieu equation for a time-dependent wave function tells us which bands are 
stable and which are unstable, but it does not address the question of how these bands are populated with time. The 
analogy with the quantum mechanical problem of a spatially periodic potential gives us insight into how it occurs. After 
the first background oscillation the x waves are not amplified at the discrete modes k which experience destructive 
interference, corresponding to the quantization (|B1I) . After the second oscillation the non-amplified modes split, and 
this splitting further continues as the number of oscillations accumulate. Asymptotically the splitting (non-amplified) 
levels cover a gap which corresponds to the stability band. 
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